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Abstract We discuss a possible interpretation of the open charm mesons DJ (2317) , D a i (2460) and the hidden 
charm mesons X(3872), Y(3940) and y(4140) as hadron molecules. Using a phenomenological Lagrangian 
approach we review the strong and radiative decays of the D* (2317) and D s i(2460) states. The X(3872) 
is assumed to consist dominantly of molecular hadronic components with an additional small admixture of a 
charmonium configuration. Determing the radiative (fJ/ip and r yif){2s)) and strong (J/ip2ir and J/ip3n) decay 
modes we show that present experimental observation is consistent with the molecular structure assumption 
of the A"(3872). Finally we give evidence for molecular interpretations of the F(3940) and V(4140) related to 
the observed strong decay modes J/ip + uj or J/if) + <f> , respectively. 
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1 Introduction 



The complexity of hadronic mass spectra allows 
for the possibility that existing and newly observed 
hadrons can possibly be interpreted as hadronic 
molecules. Such an interpretation is possible when 
the mass of the observed meson lies close to or slightly 
below the threshold of a corresponding hadronic pair. 
In the following we focus on some of the candidates 
in the heavy meson sector which are discussed and 
considered as such hadronic bound states. 

The A (3872) is one of the peculiar new meson res- 
onances which was discovered during the last years 
and where its properties cannot be simply explained 
in the context of conventional constituent quark mod- 
els. In the molecular approaches, the earliest given in 
Refs. [2, 3], it is argued that the A" (3872) can be iden- 
tified with a weakly-bound hadronic molecule whose 
constituents are D and D* mesons. This natural in- 
terpretation is due to the fact that its mass is very 
close to the D°D*° threshold and hence is in analogy 
to the deuteron — a weakly-bound state of proton 
and neutron. 

In the open charm sector the scalar D* (2317) 
and axial £) 3l (2460) mesons could be candidates for 



a scalar DK and a axial D*K molecule because of a 
relatively small binding energy of ~ 50 MeV. These 
states were discovered and confirmed just a few years 
ago by the Collaborations BABAR [4] , CLEO [5] and 
Belle ' 6 ' . The simplest interpretation of these states 
is that they are the missing j s — 1/2 (the angular 
momentum of the s-quark) members of the cs L = 1 
multiplet. However, this standard quark model sce- 
nario is in disagreement with experimental observa- 
tion since the £>* (2317) and Z? s i(2460) states are 
narrower and their masses are lower when compared 
to theoretical expectations (see e.g. discussion in 
Rcf. [7]). 

The CDF Collaboration recently reported evi- 
dence for a narrow near-threshold structure, termed 
the Y"(4140) meson, in the J/ip<j) mass spectrum 
in exclusive B + — > J/ip(f>K + decays with the mass 
m Y(4140) =4143.0±2.9(stat)±1.2(syst) MeV and nat- 
ural width r 



YX4140) 



11.7^;g(stat)±3.7(syst) MeV 



As already stressed in ' 8 ', the new structure V(4140) 
which decays to J/tp4> just above the J/tp^ threshold 
is similar to the previously discovered V(3940) [9 ' 101 
which decays to J/tpuj near this respective threshold 
Both observed states, F(4140) and V(3940), are well 
above the threshold for open charm decays. A con- 
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ventional cc charmonium interpretation is disfavored, 
since open charm decay modes would dominate, while 
the J/i/jcj) or J/tpu decay rates are essentially negligi- 
ble ' 8 ' 11 . This could be a signal for nonconventional 
structure of these Y states. 

In the following we give a compact overview of 
our recent calculations related to the strong and 
electromagnetic decay properties of the Z)* (2317), 
L> sl (2460), X(3872), y(3940) and F(4140) mesons 
interpreted as hadron molecules. Details can be found 
in the related publications I 12-21 ' on this topic. 

2 The method 

We first briefly discuss the formalism behind the 
study of hadronic molecules. As an example we con- 
sider the D*^(2317) meson as a bound state of a 
D and a K meson. We use the current results for 
the quantum numbers of isospin, spin and parity: 
I(J P ) = 0(0+) and mass m DJo = 2.3173 GeV w . 
Our framework is based on an effective interac- 
tion Lagrangian describing the coupling between the 
£>* (2317) meson and their constituents - D and K 
mesons: 

D T (x+w K y)K(x-w D y) +H.c.(l) 

where D and K are the meson doublets and lUy = 
— ^ — is a kinematic variable. The correlation func- 
tion <£>£,* characterizes the finite size of the D* (2317) 
meson as a (DK) bound state and depends on the rel- 
ative Jacobi coordinate y with x being the center of 
mass (CM) coordinate. The Fourier transform of the 
correlation function is $> D * o (p 2 E ) = exp(— p 2 E /K 2 D » ), 
where p E is the Euclidean Jacobi momentum. Here 
A D . is a size parameter, which parametrizes the dis- 
tribution of D and K mesons inside the D* molecule. 
The D* DK coupling constant g D * is determined 
by the compositeness condition ' 22 ' 23 ', which implies 
that the renormalization constant of the hadron wave 
function is set equal to zero: Z D * o = 1— ^'d* ) = 

0, where ^'d* is the derivative of the -D* meson mass 
operator. The effective Lagrangian is the basis for 
the study of the decay properties of the hadronic 
molecules. It defines the coupling of the molecule 
to its constituents but also supplies the hadron loop 
connected to the final state particles considered. The 
interaction of the constituents with external parti- 
cles is further described by effective Lagrangians ei- 
ther adjusted by data or by theory (as for example 
by heavy hadron chiral perturbation theory ' 24 ' 25 ' 



(HHCHPT)). Further details of this procedure can 
be found in Refs. [12, 21]. 

3 Strong and radiative decays of 

L> s * (2317) and L> s i(2460) 

For the Z?* (2317) we consider a possible inter- 
pretation as a hadronic molecule - a bound state of 
D and K mesons. Using an effective Lagrangian ap- 
proach we calculated the strong D* — > D s ir° and ra- 
diative D* Q — > £>*7 decays ' 14 '. A new impact related 
to the DK molecular structure of the D* (2317) me- 
son is that the presence of u(d) quarks in the D and K 
mesons gives rise to a direct strong isospin- violating 
transition D* — > D s ir° (see Fig.l) in addition to the 
decay mechanism induced by tj—tt mixing considered 
previously (see Fig. 2). 




D° D+ 

(c) (d) 

Fig. 1. Direct isospin-violating transition for 
D* s0 -> D s n° 

We show ' 14 ' that the direct transition dominates over 
the 77— 7T° mixing transition in the D* — > D a ir° decay. 




D° D+ 



Fig. 2. rj — 7T° mixing transition for D* — > D s tt° 

The radiative transition D*q — > D* s + ^ is generated by 
the loop diagrams of Fig. 3, where the last graphs have 
to be included to guarantee full gauge invariance for 
the case of a non-local vertex. 
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Fig. 3. Radiative transition D*q — > -D* + 7 

The other strong coupling vertices are hxed by 
HHCHPT or by data, for full details see Ref. [14]. 
Similar graphs apply for the strong and radiative 
transitions of the Z? s i(2460) interpreted in the present 
framework as a D*K hadronic molecule ' 13 '. Results 
for the strong decays are 

r(D* s0 ^D 3 TT°) = 47-112 keV, 
T(D sl ^Dy°) = 50-79 keV, (2) 



|X(3872)) = cos<9 

+ sin cc . 

Here 9 is the mixing angle between the hadronic 
and the charmonium components: cos 2 6 and sin 2 6 
represent the probabilities to find a hadronic and 
charmonium configuration, respectively, for the nor- 
malization 

Z D o D ,a + Z D ± D ,^ +Z Ji)U + Zj^ p = 1 . (6) 

For comparison with data we employ values for 
Z D o D ,o, Z D ± D , T , Zj u and Z Jjpp as derived in a po- 



while for the radiative modes we have 

r(D* 30 ^D* s j) = 0.47-0.63 keV, 
T(D sl ^D s -f) = 2.7-3.7 keV. (3) 

The variation in results reflects the uncertainty in the 
vertex function. The corresponding ratios of rates 
with 

R Dio = T(D* s0 ^D:y)/T(D: o ^D a ir)=OM 
R Dsl = r(D sl ^D sl )/T(D sl ^Dy) = 0.05('i) 

satisfy qualitatively the present experimental results 
of R D * < 0.059 and R D , = 0.44±0.09 Hence from 
the present observation of the radiative and strong de- 
cays of the D* and D al an interpretation as DK and 
D*K hadron molecules seems feasible. 

In Ref. [15] we further investigated the weak de- 
cays B — > +D* s0 (D sl ) , where full consistency with 
present data is achieved. This further strengthens 
the presented molecular picture. We also give predic- 
tions ' 21 ' for the weak decays involving / (980), which 
serves as a further indicator for this structure inter- 
pretation. 

4 Decay analysis of X(3872) 

The X(3872) with quantum numbers J PC = 1++ 
is considered as a composite state containing both 
molecular hadronic and a cc component. We re- 
cently showed ' 18 ' that slight binding of the DD* sys- 
tem can be achieved in a full meson-exchange model. 
Following Refs. [17, 26] we consider the X(3872) 
state as a superposition of the dominant molecular 
D°D*° component, of other hadronic configurations- 
£)±£)*=f ; J/^o;, J/i/jp, as well as of the cc charmonium 
configuration as 



(5) 



tential model. For example, for a binding energy 
e D o D ,om D o +m D *o ~m x = 0.3 MeV the estimate of 
Ref. [26] provides 

Z D o D *o — 0.92, Z D ± D ,=f — 0.033 , 
Zj^ = 0.041, Z^ p = 0.006. (7) 

The coupling of the X to its constituents, once the 
amplitudes of Eq. (7) are chosen, are determined 



7 l/2 



V2 



(\D°D* ) + \D*°D Q ))- 



yl/2 



V2 



(\D+D*-) + \D-D*+)) + Zy;jJ^)+zX 2 p \J^ P ) 
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again by the compositeness condition ' 16, 17, 19 '. 

A new measurement by the BABAR Collabora- 
tion gives clear evidence for a strong radiative decay 
mode involving the ip(2S) ' 27 ' . They indicate the mea- 
sured ratio of 



B(X(3872)^^(2S)j) 
£(A(3872)^ J/i/yy) 



= 3.5±1.4. 



(8) 



As known from previous calculations in the molecular 
approach the radiative decay width of the A (3872) — > 
^(25)7 is always smaller than the one involving 
J/-07- It is therefore expected that the ratio of Eq. (8) 
gives some constraint on a possible charmonium com- 
ponent in the A (3872). 

In Fig. 4 we display the diagrams relevant for the 
radiative decays A — > J/ip^f and A — > i/j(2S)j. The 
last diagram of Fig. 4 which is generated by the J/ipV 
(V = p,uj) component of the A(3872) only contributes 
to the J/ipJ final state. 



60 keV and is much larger than the one for the ip(2S)j 
channel (about 0.3 keV). Inclusion of a charmonium 
component in the X(3872) configuration leads to the 
results of Fig. 5 indicating the radiative decay widths 
in dependence on the cc admixture. 




-sin9 



Fig. 5. Radiative decay widths of the X(3872). 
Solid and dotted curves are the results for the 
transitions X — > J/tp"/ and X — > tp(2S)^/, re- 
spectively. 





Fig. 4. Diagrams contributing to the radiative 
transitions X(3872) -> J/ip+-y and X(3872) -> 
V>(2S)+ 7 . 

Results ' 16 ' for the radiative decay modes involving 
the hadronic components of the A(3872) are con- 
tained in Table 1. The value of the binding energy 
is set to 0.3 MeV with the probabilities taken from 
Eq. (7). 

Table 1. Radiative decays of hadronic components. 



Hadronic config. 


r(X(3872) ->7J/i/>, 7^(25)) keV 


DD* 
J/ipV 

DD* and J/ipV 


60 - 12O(J/i/0 0.3 0(250) 
6(J/V) (V(2S)) 
30 - 65 (J/V>) 0.3 (ip(2S)) 



The inclusion of the J/ipV components in the 
structure of the A (3872) leads to a decrease in the 
A — > J/V>7 rate because ip(2s)V components are ab- 
sent in the X(3872). Note that taking into account 
the hadronic components only the ratio of Eq. (8) 
cannot be explained. If no cc component is present 
in the A (3872) the decay width for J/V>7 is about 



A small admixture of a cc component is essential 
to reproduce the large measured value for R 2 s ' 16 '- 
This feature is due to the destructive interference be- 
tween the J / V>7 decay amplitudes arising from the cc 
and the hadronic components. For mixing values of 
about sin# rj —0.16 in the case of e = 0.3 MeV the 
decay width for ijj{2S)^ may exceed the one of J/ip'y 
and a consistent explanation for the observed ratio 
i?2s is obtained. 

To get estimates ' 17 ' for the decay widths of 
A(3872) -» J/i/j + h with h = 7r+7r-7r ,7r + 7rV°7,7 
we use the results of Ref. [28], which are based on 
the assumption that these decays proceed through 
the processes A to J/tpui and J/ipp- Here we esti- 
mate ' 17 ' both short and long-distances effects only for 
the A — > jJ/'ip decays using our previous result, while 
the A — > J/ip+h decays only take into account short- 
distance effects. Results for the J/ip decay modes are 
given in Table 2. 

Table 2. Properties of X — > J/ip + h decays. 
Numbers in brackets are valid for explicit val- 
ues for Zj^p and Zj. u of Eq. (7). Data for 
ratios are from Refs. [27, 29]. 



Quantity 


Nonlocal case 


r(X-j. J/lpTT+TT-), keV 


9.0xl0 3 Z^ p (54.0) 


F(X-> J/ipn+TT-iT ), keV 


1.38xl0 3 Z.^„ (56.6) 


F(X-> J/i>7T -y), keV 


0.23 x 10 3 Zj^ (9.4) 


r(Jf^J/i/)7r + -/r-7T U ) 
r(X->.J/i/.ir + 7r-) 

= 1.0±0.4±0.3 


1.05 


r(Jf-i-J/i/>7r + 7r-) 

=0.14±0.05; 0.33±0.12 


0.10 
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Present data for the ratio of rates ^tttt + > 
indicate a strong isospin violating J/tpTr + TT~TT° decay 
mode. Similarly, the measured ratio r (x^j/i!J+-l-) 
points to a large radiative decay channel. Both ra- 
tios of rates can be fully reproduced in the molec- 
ular picture. Note that all the decay channels con- 
sidered so far are fed by the subleading J/i/jlj, J/tpp 
and cc components. To work out the influence of 
the leading molecular components we consider the 
hadronic X — > Xu + (tt° , 27r) and X — > J/iJj + (2ir, 3ir) 
decays ' 17 ' . Here the values of J = 0, 1, 2 correspond to 
the J p = + , 1 + , 2 + quantum numbers of the charmo- 
nium states. Because of the dominance of the D°D*° 
component in the transitions of X into charmonium 
states Xcj and pions we estimate these decays using 
only this component. The two-body decays are de- 
scribed by the (D°D*°) loop diagram shown in Fig. 6. 




Fig. 6. Diagrams contributing to the hadronic 
transitions X(3872) -»XcJ + ""°- 



Inclusion of the charged (D ± D* T ) loops gives a 
further correction to the decay widths. The possible 
couplings of D(D*) mesons to pions and charmonia 
states is taken from HHCHPT. Similar graphs ap- 
ply for the three-body decay modes. Results for the 
hadronic decays involving the neutral (D°D*°) com- 
ponent only and the full results are indicated in Ta- 
ble 3. 



Table 3. Properties of X 
Here Zo = Z D o D *o. 



XcJ + niv decays. 



Quantity 


D°D*° loop 


D°D*° + 
+D-D*+ 


r(x^ X co + TT°), kev 


41.1 Zo (37.8) 


61.0 


r(X-> X cO + 27r°), eV 


63.3 Z D (58.2) 


94.0 


r(X^ X ci + 7r°), keV 


11.1 Z D (10.2) 


16.4 


r(X->Xci+27r°), eV 


743 Zo (683.6) 


1095.2 


r(X^Xc2 + 7r°), keV 


15 Z D (13.8) 


22.1 


r(X->Xc2 + 27r°), eV 


20.6 Zo (19.0) 


30.4 



Note that explicit numbers refer to a binding en- 
ergy of 0.3 MeV with the probabilities taken from Eq. 
(7). The decay pattern of Table 3 involving pions 
and Xcj states is sensitive to the leading molecular 
(D°D*°) component. These predictions ' 17 ' can serve 
to possibly identify the full hadronic composition of 
the X(3872) in running and planned experiments. 



J/W(= 

Y(4140) 



to, < 



decays of Y~(3940) and 



It was suggested in Ref. [30] that both the Y (3940) 
and Y(4140) are hadronic molecules. These hadron 
bound states can have quantum numbers J PC = 
++ or 2 ++ whose constituents are the vector charm 
D*(D*) mesons: 



|Y(3940)) = 
|y(4140)) = 



1 

\d:+d 



D*+D*-) + \D*°D* )) 



(9) 



For the observed F(3940) and F(4140) states we 
adopt the convention that the spin and parity quan- 
tum numbers of both states are J PC = ++ . The 
coupling of the scalar molecular states to their con- 
stituents is again set up by the compositeness condi- 
tion. To determine the strong Y — > J/tpV and two- 
photon Y — >• 77 decays we have to include the cou- 
plings of D*(D* g ) mesons to vector mesons (J/ip, w, 
4>) and to photons. The couplings of J/tp, to, <j> to 
vector D*(D* a ) mesons are taken from the HHChPT 
Lagrangian ' 24, 25 ' . The leading-order process rele- 
vant for the strong decays Y(3940) — > J/ipu and 
Y"(4140) — > J/ip<fi is the diagram of Fig. 7 involving 
the vector mesons D* or D* in the loop. 




Fig. 7. Leading order diagram for Y — s- J/ipV 
decay. 

We also consider the radiative Y"(3940) /Y(4140) -> 
77 decay widths ' 20 ', which serve as a prediction for 
possibly identifying the molecular structure of these 
Y states. 




Fig. 8. Diagrams relevant for the radiative de- 
cay Y — > 77. 
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The relevant diagrams which arise from coupling of 
the charged D* ± (D* ± ) mesons to photons and from 
gauging the strong interaction Lagrangian are dis- 
played in Fig. 8. The numerical results for the quan- 
tities characterizing the strong J/ipV (V — lj, (f>) and 
radiative 27 decays of Y(3940) and Y(4140) are con- 
tained in Table 4. 



Table 4. Decay properties of Y(3940) and 
Y(4140). 



Quantity 


Y(3940) 


Y(4140) 


r(y-> j/ipv), Mev 


5.47 ± 0.34 


3.26 ± 0.21 


r(Y-»77), kcV 


0.33 ± 0.01 


0.63 ± 0.01 


r(y -»■ j/w) 


0.61±0.06 


1.93±0.16 



The predictions of r(V(3940) J/ipLu) = 

5.47 MeV and r(Y(4140) -> J/tp(p) = 3.26 MeV for 
the observed decay modes are sizable and fully con- 
sistent with the upper limits set by present data on 
the total widths. The result for r(Y(3940) -> J/ipu) 
is also consistent with the lower limit of about 1 
MeV ' n '. Values of a few MeV for these decay widths 
naturally arise in the hadronic molecule interpreta- 
tion of the Y(3940) and Y"(4140), whereas in a con- 
ventional charmonium interpretation the J/tpV de- 
cays are strongly suppressed by the Okubo, Zweig 



and Iizuka rule ' . Further tests of the presented 
scenario concern the two-photon decay widths, which 
we predict to be of the order of 1 keV ' 20 ' . Results 
for the strong J/ip decays are quite similar for the 
J PC = 2++ assignment, which cannot be ruled out at 
this point. 

6 Conclusions 

The approach to hadron molecules based on the 
compositeness condition constitutes a consistent field 
theoretical tool. The determination of the decay 
properties of the open charm mesons D* (2317), 
D s i(2460) and the hidden charm mesons X(3872), 
Y(3940) and Y(4140) in comparison to present data 
is fully in line with a molecular interpretation of these 
states. We also give further predictions of decay prop- 
erties which can be tested in future experiments. Cur- 
rent results are very encouraging in ultimately iden- 
tifying hadronic molecules in the meson spectrum. 

We thank Yubing Dong, Yong-Liang Ma and 
Sergey Kovalenko for extensive collaborations on the 
topic of hadron molecules. V.E.L. is on leave from 
the Department of Physics, Tomsk State University, 
634-050 Tomsk, Russia. 
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